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Abstract 



A one-dimensional system of nonintersecting Brownian particles is constructed as the 
diffusion scaling limit of Fisher's vicious random walk model. N Brownian particles start 
from the origin at time t = and undergo mutually avoiding motion until a finite time t = 
T. Dynamical correlation functions among the walkers are exactly evaluated in the case 
with a wall at the origin. Taking an asymptotic limit iV — > oo, we observe discontinuous 
transitions in the dynamical correlations. It is further shown that the vicious walk model 
with a wall is equivalent to a parametric random matrix model describing the crossover 
between the Bogoliubov-deGennes universality classes. 
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1 Introduction 



The vicious walk model was first introduced by Fisher and applied to wetting and melt- 
ing phenomena P3 12]. It recently attracts renewed interest in statistical and mathematical 
physics, since intimate connections were established to other research fields, such as Young 
tableaux in combinatorics [31 HJ Ej, asymmetric exclusion process (ASEP) in nonequilib- 
rium statistical mechanics jH], Kardar-Parisi-Zhang (KPZ) universality in surface growth 
pro cess [7JIHJII3 an d the theory of random matrices |1U| II lj. In the context of random matrix 
theory, the ensembles of vicious walkers in one dimension correspond to discretizations of 
random matrix ensembles. 

Suppose that there are iV walkers on a lattice Z = {• • • , —2, —1, 0, 1, 2, • • •}. Starting 
from N distinct (even integer) positions 2s\ < 2s 2 < ■ ■ ■ < 2s^, at each time step each 
walker moves to the left or right one lattice site with equal probability. Let us denote 
the position of the j-th walker at time k > by R s k 3 . Walkers are "vicious" so that two 
or more walkers are prohibited to arrive at the same site simultaneously. That is, the 
nonintersecting condition 

R£ < R s k 2 < ■ ■ ■ < R s k N , l<Vk<K (1.1) 

is imposed. In this paper, we further impose a condition 

Rk > 0, 1 < VA; < K, (1.2) 

which implies that there is a wall at the origin. Let us define V{R S ^ = 2dj) as the 
realization probability that the vicious walkers arrive at (even integer) positions 2d\ < 
2d,2 < ■ ■ ■ < 2d,N at an even integer time K. Utilizing the Lindstrom-Gessel-Viennot 
theorem and the reflection principle of random walks, we find an explicit formula [Sj 



V(B% = 2dj) = 2~ KN det 



K \ ( K 



K/2 + s k - djj \K/2 + s k + dj + 1 / 



• (1-3) 

j,k=l,2,-,N 



This realization probability can be simplified in the diffusion scaling limit |12t ITSl ITU Ho] . 

Let us introduce a positive number L and set K = Lt, Sj = y/Zxj/2 and dj = \fLyj/2. 
Then we can easily find 



Vi, V2,---,Vn\ xi%2, ■■■ ,x N ) = lim ( ^f) V(R^ X ° /2 = ^/Tyj) 

L^oo \ 2 J 

(2nt)~ N/2 det exp l~-^( x j ~ ^fe) 2 } - ex P [~^ x i + 2/ fc ) 2 } 



This function gives the nonintersecting probability of the Brownian particles on the 
rescaled lattice Z/(V^/2) up to times t depending on the initial positions {xj} and final 
positions {%•}. Therefore the probability amplitude that the vicious walkers are located 
at x{, xl, ■ ■ ■ , x 3 N at times t,-, j — 0, 1, 2, ■ • • , M, is given by 

P{to'-> x \t x 2i ' ' ' i X N) ^1) X li X 2i ' ' ' i %i ' ' ' ) ^M] %l i%2 i'''i x n) 

AI-l 

= II ^ (tj'i x i' x 25 ' ' ' 5 X 7V) x i j x 2 r - ') 2 ^ )> (1-5) 

3=0 
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where 



ip (s; xi, x 2 , ■ ■ ■ , x N ; t; y u y 2 ,---, Vn) 
= f(t - s; yi, y 2 , ■ ■ ■ , Vn | x u x 2 , ■ ■ ■ , x N ) 



A/"(T -t;y 1 ,y 2 ,---,y N ) 
Af(T - s;xi,x 2 , ■ ■ ■ ,x N )' 



(1.6) 



A/"(t; xi, x 2 , • • • , x N ) = / dyidy 2 • • • dy N f{t; y x , y 2 , ■ ■ ■ , y N | x u x 2 , ■ ■ ■ , x N ). 

J0<yi<y2<—<yN«x> 

(1.7) 

The dynamical correlation functions among the walkers at times ti, t 2 , ■ • • , tM are defined 

as 

( 1 ... 1 2 ... 2 . . . . M ... M \ 
= T\dx°i TT drrj • • • / TT dx 



X 



j=ni+l 

Po({*?}) II ^ T (t m ,{^};t m+1 ,{a;f +1 }) 



M-l 



1.8) 



m=0 



Here Po({^j}) is the initial distribution at t = 0. Let us suppose that all the Brownian 
particles start at the origin so that we can set Po({x^}) = Hj S(x°). Making replacements 



x ) ~^ yy}' x ] ~^ \fy]> " * ' x f 



1.9) 



we rewrite the dynamical correlation functions as 

P(Vli ' ' ' j ; ' ' ' j 2/n 2 ' ' ' ' ' fl ' ' ' ' ' ^n M ) 



oc 



roc /-oo fOO foo 

/ n dyj/ n dyj---/ n 

JU i=ni+l JU 7=n 2 +l J0 i=n M +l JU .7=1 



M+l 
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j=m+l j=n 2 +l j=n M +l 

N N M 

x mv)- vi) n s ^f +1 - ^ M+1 ) n det kw, 

j>fc j>fe m=l 



j,k=l,2,-,N 



(1.10) 



Here 



^(^y) 



x 



2 v / 2^T ^27r(t 2 - ti)y 
1 



cxp 



{"277"} 



exp 



2(i 2 -ti) 



(Vx - v/y) 2 [ - exp 



2(i 2 -ti) 



(1.11) 



and 



x 



9 m (x,y) 



exp 



2^2n(t m+1 - t m )y 
1 



2(iyn+l ^r, 

2 < m < M. 



(-v/x - ^fyf \ - exp 



1 



■(v^ + W) 2 
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In the case with no wall, the vicious walk model in the diffusion scaling limit is equiva- 
lent to the eigenvalue dynamics of parametric random matrices belonging to the standard 
symmetry class |13| Hlj . Similarly, in the presence of a wall, it will be shown that the 
vicious walk model and parametric random matrices with the Bogoliubov-deGennes sym- 
metry are equivalent. The Bogoliubov-deGennes matrix model was proposed by Altland 
and Zirnbauer as an effective model of mesoscopic normalconducting-superconducting 
hybrid structures [TBI El- 

This paper is organized as follows. In §2, quaternion determinant expressions for the 
dynamical correlation functions are presented. In §3, we rewrite the quaternion determi- 
nant expressions in terms of the Laguerre polynomials. In §4, asymptotic forms of the 
dynamical correlation functions are evaluated in the limit iV — > oo. In §5, an equivalence 
between the vicious walk and the Bogoliubov-deGennes matrix model is demonstrated. 



2 Dynamical Correlation Functions 

2.1 Quaternion Determinant Expressions 

We begin with the definition of a quaternion determinant [T5| . A quaternion is defined as 
a linear combination of four basic units {1, ei, e 2 , 63}: 

q = go + qiei + g 2 e 2 + q^. (2.1) 

Here qo,qi,q2 and qs are real or complex numbers. We call go the scalar part of q. The 
quaternion multiplication is associative but in general not commutative. The multiplica- 
tion rule of the four basic units are given by 

1-1 = 1, 1 • ej = ej • 1 = ej, j = 1, 2, 3, (2.2) 

e i = e l = e l = eie 2 e 3 = -1. (2.3) 
We define a dual q a quaternion q 

q = qo - qiei - g 2 e 2 - <?3e 3 - (2.4) 

A dual is an analogue of the complex conjugate of a complex number. For a matrix Q 
with quaternion elements qji, we can also define a dual matrix Q = [qij]. The quaternion 
units can be represented as 2 x 2 matrices 

-1 " 

1 ' 





-i " 




" % 


e 2 -> 


-i 


, e 3 -> 


-i 



(2.5) 



Let us now introduce a quaternion determinant Tdet. For a self-dual Q ( i.e., Q = Q 
), it is defined as 

i 

Tdet Q = E(" 1 )' V " i U(qa b q b c • ■ ■ q da )o- (2.6) 
p l 



1 
1 



ei 
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Here P denotes any permutation of the indices (1, 2, • • • , N) consisting of I exclusive cycles 
of the form (a — > 6 — ► c d — > a) and (— 1) is the parity of P. The subscript 

means that the scalar part of the product is taken over each cycle. If all the elements of Q 
have only the scalar parts, every element is commutable so that a quaternion determinant 
becomes an ordinary determinant. 



In both cases N even and odd, we define quaternion matrices B^ , //, u — 1, 2, 
in terms of the following 2x2 representations of the quaternion elements: 



M 



jl J I 



J, I 



1,2, 



N. 



(2.7) 



The matrices S^ U ,I^ U and are defined in the following subsections. Applying the 
integration technique developed in Ref-jTHj to the integral (jl.K)j) . we find the quaternion 
determinant expression of the dynamical correlation functions 



p(y{, 

(i,u-- 



> Vni ' Vl ' 



> Unit 



■v M 



v M N 

' t>n M ' 



Tdet[S' B/ (n M ,n w )], 



1,2, 



M, 



(2.8) 



where each block B^ v \ 
and n v + 1, n v + 2, • • • 



is obtained by removing the n M + 1, n M + 2, 
iV-th columns from B^ u . 



iV-th rows 



2.2 Skew Orthogonal Polynomials 



In order to define the matrices S^ u , and Dj t , we need to first introduce skew orthog- 
onal polynomials. In terms of 

G mn (x,y) 

5(x — y), m = n, 
g m (x, y), m = n — 1, 

/ dy m+1 dy m+2 • • • d?/ n _i 
Jo 

X5 ,m (^,l/m+l)fl' m+1 (l/m+l ) ym+2) ' ' ■ (y n -U y) , m < U - 1, 



we define 



(2.9) 



F mn {x,y) = / dz' dz{G m M+1 {x,z)G n M+1 (y,z') 
Jo Jo 

- G n M+1 (y,z)G mM+1 (x,z')}. (2.10) 



In terms of an antisymmetric inner product 

(f(x),g(y)) m 



1 POO POO 

--- dx dyF mm (y,x)[f(y)g(x)-f(x)g(y)}, (2.11) 

2 jo Jo 

monic polynomials R\(x) = x k + ■ ■ ■ of degree k are constructed so that they satisfy the 
skew orthogonality relations: 

(</)>! = -(R 1 2l+1 (x),R 1 2j (y)) 1 = rj 5 jh 
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(i2j,-(s),i&(y))i = 0, {Bk j+1 {x),R\ l+1 {y))i = 0. (2.12) 
Defining a set of functions R™(x), m = 2,3,---,M + las 

RT(x)= / d^(y)G lm (|/,x), (2.13) 

we can immediately find the following skew orthogonality relations for m — 1, 2, • • • , M+l: 

<i?™.(x),i^ +1 (y)) m = -<i^ +1 (x),i^.(y)> m = r^ 7 , 

(i^-(x),i^(y)> m = 0, (i^- + i(x),i^ +1 (y)> m = 0. (2.14) 
At this stage we need to consider the cases N even and odd separately. 

2.3 The Case JV Even 

Let us now introduce matrices D mn , I mn , S mn , F mn and G mn as 

(N/2)-l 1 

^r= E -K(yr)^ + i(z/r)-^ + i(yr)^(yr)], (2.15) 

fc=0 ' k 
(JV/2J-1 1 

T = - E -K(yr)^ + i(yD-^ + i(yr)^(yD], (2.16) 

fc=0 rfc 

(AT/2)-l j 

?r = E -K(yr)^ + i(z/r)-^ + i(yr)^(z/r)], (2.17) 



k=0 



where 

Further we define 



and 



PCX) 

*?(x)= / F™(y,x)R™(y)dy. (2.18) 

FT = Fmn (v?>vi)> ( 2 - 19 ) 

(™» = 0, m>n, 

G rnn = ^rn ^ m<n {22Q) 



jmn jmn _|_ ptnn Qimn grnn Qnya fc^ € 2^\\ 

3^ 3^ ' 3^ 3^ 3^ ^ ' 



2.4 The Case iV Odd 

In terms of skew orthogonal functions R™(x), another set of functions are defined as 

i%(x)=R?(x)--?2-B%_ 1 (x), n = 0,1,.-., N -2, (2.22) 

AjUC*) = i^s), (2.23) 

where _ 

<r fc = / dxf M (x)R»? (x) (2.24) 
Jo 



with 



f m (x)= / dyG m0 (x,y). 
Jo 



Using the definitions 



f? = riv?) 



and 



1 



we introduce the matrices D mn , I mn and S mn in the case N odd as 

(JV-3)/2 1 

e -K(yr)% + i(yD - ^(vrAwi. 



Y)mn 



k=0 Tk 

(N-3)/2 



fc=0 



and 



where 



+ 



(7V-3)/2 



CTN-1 



cjv-i 



s 1 : 



1 



jl 



E -[^(2/D^Vi(% n ) - *%+i(v?)Rh{v?)] + fT R > 

k=0 Tk 



n 
I ! 



(y,x)R%(y)dy. 



Then the matrices I mn and S' mn in the case N odd are given by 

fmn jmn i rpmn Qmn omn /~imn 

'jl ~ 'jl + '' ;/ ' ^jl ~ °jl ( 'J< ' 

where matrices F mn and G mn are defined in eqs. (|2.19|) and (|2.2(J|) . 



(2.25) 

(2.26) 
(2.27) 

(2.28) 



(2.29) 

(2.30) 
(2.31) 

(2.32) 



3 Description in terms of the Laguerre Polynomials 

In order to derive the asymptotic behavior of the dynamical correlation functions, it is 
convenient to rewrite the quaternion determinant formula presented in §2 in terms of the 
Laguerre polynomials. Note that (£m+i = T) 

G ln (x,y) 

i i r i 



2V2rtiJ 2 7r(t n -t 1 )y eXV \ 2*1* 



x 



exp 



2 (i n — ii 
1 1 1 



(Vx - jyf \ - exp 



2(tn - h) 



v / 2^ii'^2 7 r(i n -i 1 ) 2 v / ^ 
1/1 1 



x exp 



exp 



2 Vii i n — ii 
1 



x > exp 



2(i„ — ii) ' 



in ii 



xy > — exp 



i« — ii 
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(3.1) 



Let us define 

and rescale x and y as 

G ln ( Cl X,c n X') 



tni^T t r . 

T 



(3.2) 



1 



y/^tij2n(t. 



^exp ( — (A + X')\ exp |%A' 



x 



exp 
1 



2 V / ^AA 7 
1 - Z n . 



— exp 



2 V ^AA 7 ' 
1 - z n 



\j^T\j Z n{^- Zn) 



x A 1/2 e" A exp 



2T 



A' 



unj^m Ll (A)Lj (A) - 



E 



Here 



Zn. 



t\ IT — t n 



(3.3) 



(3.4) 



IT — ti t n 

and L^ 2 \x) are the Laguerre polynomials (with a parameter 1/2). For m > 1, we can 
similarly derive 



G mn (c m X, c n X') 



1 27r(t m - ti) z n (l - Zn) 



y/C m C n \| 27T(t„ - ti) \| z m (l 

Zm ) 



x A 1/2 exp 



A 

~2T 



exp 



1 - ^ ) A' 
2T 



r (j + ( 3/2))^ (A) ^' (A) 



I" 

j=0 



Let us now define monic orthogonal polynomials 

C n (x) = {-l) n n\L^ 2 \x) = x n + ■ ■ 

satisfying 

/ x l ^ 2 e- x C j (x)C l (x)dx = hjSji 

J 

with 

^ = (j + (1/2))!j!. 
We then introduce skew orthogonal polynomials 



where 



.(2ife)! (2fc-r- (1/2))! 



(3.5) 

(3.6) 
(3.7) 
(3.8) 

(3.9) 



(-1)' 



v^F (2Jfe -r)!H 



«2fe+l r — ( — 1) 



.(2k + 1)! 



7T 



2A;-r + 



1\ (2fc-r- (3/2))! 



4/ (2fc-r + l)!r! 



(3.10) 



and 



Xr. 



Zm+1 



It is known that BJ^{x) satisfy skew orthogonality relations |2U1 |2~T] 

((R%(x),R™ l+1 (y))) m = -((R™ l+1 (x),R%(y))) m = ff5 jh 
((R%(x),R%(y))) m = 0, ((Rf j+1 (^Ki+i(y)))m = 



with 



r7 = 4(2j + l)!(2j)!(x 
Here the antisymmetric inner product is defined as 



«/(*), 0(y)», 



dx r dyx^Y^e-^e-y^F^iy, x) 
Jo Jo 



= o / dx \ 
2 Jo Jo 

x [f(y)g(x) - f(x)g(y)} 



where 



F mn (x,y) 



o Jo 



dz' / dz- 



zz') 1 ^ 



with 



x [g m M +\x, z)g n M+1 (y, z') - «f M+ \y, z)g m M+1 (x, z')} 
g mn (x,y) = x^y^e-^e-y^Y, 



, z. T . 



j-0 h 3 



We can readily find (m > 1) 

G 1M+1 (c 1 A,c M+ iA / ) = 
G mM+1 (c m X,c M+1 \') 



A 



1/4 



V2^T2T-t 1 \X' 
1 



— e 



- A/ V M+1 (A,A'), 



t m t\ 



\fC m CM+l V t-M+1 ~ t\\ z miX ~ z rn) 



X 



A 



1/4 



exp 



1 - 



T / 2 



3 mM+1 (A,A') 



so that 



F n (c 1 \,c 1 X') = 
F mm (c m \, c m \') 



1 / ^y(AA , ) 1/4 e- A /V A '/ 2 J p 11 (A,A'), 



x (AA') 1/4 exp 



1 - 



2ttT Vci 

C*A/+1 — ^1 ZM+l{l — ZM+l) 

Cm ^M+l tl Z m {\ Zfn) 

t<m \ A 



T 2 



exp 



U 7-7") \ /\ 



T / 2 



F mm (A,A') 



It is now straightforward to see that the polynomials 

Rl(x) = (c^RKx/c,) 
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satisfy the skew orthogonality relation ()2.12|) with 

: (2j + l)!(2j)! 



2._ „...({t x yV j+2 



71 



T 



since 



((R?(x),Rr(y)}} m = 2Tr^(c 1 z n 



-^ l {Rf(x),Rr(y))r 



(3.20) 



(3-21) 



Then, from the definition, we obtain (m > 1) 

1 



R?(x) = {z m f{ Cl ) 
t, 



k+1 



7T 



x exp 



1 



277 a 



x 



RT(x/c m ). 



^"m ( 1 Z m ) 



(3.22) 



Putting the above results into the quaternion determinant formula yields 

' ) Vn\ 1 Vl i ' ' ' 1 Vri2 1' ' ' iVl 5 ' ' ' ) Vn M ' 

= Tdet [B^(n^n u )] 

M 

= n(Q)^Tdet[^(n M ,^)], 
i=i 



(3.23) 



where each block 5 Ml/ (n M , n v ) is obtained by removing the n M +l, n^ + 2, • • • , iV-th rows and 
n v + 1, n v + 2, • • • , iV-th columns from B^ v . The quaternion elements 5^ are represented 
in terms of the scaled variables 

K = y?/ c m (3-24) 



as 



j,l = l,2,---,N. 



(3.25) 



The definitions of S fll/ {x 1 y) ) I^ix, y) and D^ v {x,y) are (for even N) 



(N/2)-l 

E 



k=0 1 k 



(XrnY k R? k (x)( X n) 2k+1 R n 2k+1 (y) 



(xJ 2k+1 R2 k+ i(x)(xn) 2k R n 2k (y) 



(3.26) 



{x,y) 



(JV/2)-l 



E 

fc=0 



r 



2fc+l 



(y) 



- (Xr, 



\2fc+l ( f 1 m 



^+i(^)(xn) 2fc+1 ^(y) +F mn (x lV ), 



(3.27) 



5 m "(x,y) 



(JV/2)-l 



k=0 



~M+1 



(x m ) 2k ^(x)( Xn r +L R n 2k+1 (y) 



2k + l T3TL 



{Xm) 



2^+1 <^> m 



x){Xn? k R n 2k {y) ~G mn (x,y), 



2fc+l 



(3.28) 
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where _ 

- r°° - I,- 

<f>Z{x)= F mm (y,x)y 1 / 4 e-y/ 2 R^(y)dy (3.29) 
Jo 

and 

G mn {x,y) = 0, m>n, 

G mn (x,y) = g mn {x,y), m < n. (3.30) 



4 Asymptotic Correlations 

Let us consider the asymptotic limit N — > oo of the dynamical correlation functions. A 
new result should be searched in the neighborhood of the origin, since, in the regions 
far from the origin, the asymptotic behavior of the dynamical correlations should not be 
changed by the presence of the wall. In order to see the asymptotic correlations around 
the origin, we define scaled temporal and spatial variables v m and X™ as 



y m X m 
c m > ' ^ N 



Taking the expansion ()3.9|) and putting an asymptotic formula 



lim 4 a) (x) ~ e x/2 (-) a/2 J a {2^x~) (4.2) 



n^oo 



.X, 

(J a (x) is the Bessel function) which holds uniformly for x = 0(l/n), we can readily derive 

(\ ) 2k AT/9 3 / 4 fi tj 1 / 4 / 

*5<*W W ~ V^l d "(T^^(2V^)e»-, (4.3) 

where 

e = — . (4.4) 

The skew orthogonal polynomials with odd order are rewritten as 

IT (x) {Xm)2k+1 - ( 2k + W 2 W L M)( x ) (2fc-(i/2))! (1/2) 

- 7E (2 1: ' + | 1 (' ))! [4 1/2) (^)(^)- - L^\ X ){ Xm r\ (4.5) 
^2 (2A;-r + l)! 

Here the contribution to the asymptotic form comes from the last term and is evaluated 

as 



1 \2fc+l o r l 1 (I i 



(2fc + i)i ^(xey/^Jo '(i-T]y/ 2 dT] 



■ (4-6) 



11 



We substitute the asymptotic forms of R™(x) into ()3.26|) to derive 

/V 3 / 2 rl 

D mn (x/N,y/N) / d^ 1/2 

47T JO 

/•i r ? 1 / 4 J 1/2 (2 v ^)e^" /•! 1 d f 1/4 /— ^ 

- 7o dr? (T^F 2 " 7o de (wF^ ( e J ^(W)e 

We now introduce the inverse matrix [fiji\j,i=i,-,k of [c*!jz]j > i=i } —,& as 

c fc (x) = (x m )~ fc E^r(^)(x m ) r , 



r=0 



where 

Pit 2r 



(fc-2r- (3/2))! fc! 
2^ (fc - 2r)! (27)! 



(-l) fc fc! [0-2r-i)/2] ( fc _ 2r - 2Z - (5/2))! 
^ 2r+1 ~ 2^ (2r + l)! ^ (fc-2r-2Z-l)! " 

In terms of the inverse matrix /3ji, the function $>™(x) is written as 



00 



^(^(X m ) 2fc = x^e-^rr 1 £ 2fc+1 . 

i/=2fc+l 



Using the above expression and identities 



2/--^! 00 (21-- 



,=o < 2 ' + !> ! i-lr+i (2' + 1)! ' 

f( 2 '-D ! f ( 2 '-D ! 

a (2')! ,ir+, (20! ' 

we can rewrite the function &™ k {x) as 

^ (x) w x_1/v/2 

" ^ r= t +1 (2r+(l/2))!^ ^ (2/ + 1)! 



L V (2r~l)! r (i/2 )fT w v r2r+1 v (2/ ~(3/2))! 



^i(2r-(l/2))! — /VAm/ ^ (2/)! 

12 



+ 2 Z ' 2 J2 

r=k+l 



(2r)! 



Lg /2 \x) (2r-l)! L^ 2 }(x) 



(2r + (l/2))! (xj* (2r-(l/2))!( Xm )2r-i 



The first and second terms contribute to the asymptotic limit and give 

2 ,„„,i, A f°° 1 

' I T] 1 / 4 ^- 1)1/2 



C \2fc 



(2A;)! v^F 
Then we put $2fc + i(a;) in a form 



JiMe V x)e- ai 



(4.12) 



(4.13) 



(Xr 



\2fc+l 



X 



-1/4 x/2 



2,+lV ~ ; (2A; + 1)! 4 

(2r-2Jfe-(5/2))l (2r-2)! 



E 



(2r-2A;-2)! (2r-(l/2))! 

x (Xm)- 2r+1 [(2r - l)2#2(s) - (2r - I) L^(x) 

~ (2r-2fc-(3/2))! (2r - 1)! 
+ (2r-2fc-l)! (2r + (l/2))! 



X (Xm) 



-2r+l 



(2r)4r /2) (a:)-(2r + |)4;.j;(.r) 



1 



(1/2), 



x (2r-2Jfe-3)! (2r-(3/2))f 



-2r+l 



r=k+l 



{2r-2k-2)\ (2r-(3/2))! 



-2r+3l 



v-2r+2n 



(4.14) 



and substitute (|4.2|) to derive 

(Xrn) 2k+1 4 



2fc+l 



>/iV) 



(2fc + l)! v^F(^) 3 / 4 7i '(77 — l) 1 /2d77 



drj- 



71/4 



J 1/2 (2j6r)x)e 



-9rjv„ 



Putting the expansion ()4.8|) into ()3.15|) and using (|4.10|) yields 



(4.15) 



00 1 

E 



fe=0 r fe 



A/+1 



(X m ) 2fc $ 2 rt fc (x)(Xn) 2fe+1 ^ fc+ l(l/) 



2fc+l ( f 1 n 



(Xr. 



\2/c+l/F>m 



»(Xn) 2fc+1 ^ fc (l/) 



2fc+l 



(4.16) 



so that 



E 

k=N/2 ' k 



{x m ) 2k+1 ^ k A^Xn) 2k+l n k {y) 



(4.17) 



Asymptotic forms ()4.13|) and ()4.15|) are substituted into the above expression and yield 



(x/N,y/N) 



ATl/2 



7T Jl 



d9 
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A ^ ^f. - 1)1/2 / d fe)^d7 7I74 J V2(2^)e 



^1/4(^-1)1/2 ^ S (^- 1)1/2 [^1/4 

00 J l/2 {2^Wqy)z~ eriVn f°° 



Moreover we put (jOjl . (ELTfl) . (jHT3j) and (ETToT) into (ET^Kj) and find 

S m "'(x/iV, y/N) + G m "(x/iV, y/JV) f'dO 

271 Jo 

- r J 1/2 (2y/^)e-^ ^ 1 d J 1/4 /— ^ 

x A ^ ^1/4^-1)1/2 A de (w)^de{ e ,J wM e tv)*' 
- A d ^(^iy^^i^ Jl / 2(2 v^ )e |A de — a-o 1 / 2 



(4.18) 



(4.19) 



The function G mn (x,y) is identical to g mn (x,y) when m < n. The asymptotic limit of 
g mn (x,y) is derived from ()3.16|) and (|4.2j) as 



g mn (x/N,y/N)~N dsJ 1/2 (2^)J 1/2 (2^)e-^- v ^ s . (4.20) 



Substituting (|T7 jl . (fOB)l . (jO?)l and (TO0jl into JH23|), we can see how the dynamical 
correlation functions asymptotically depend on the scaled variables v m and X™. 

5 Bogoliubov-deGennes Matrix Model 

In this last section we show an equivalence relation between the vicious walk model with a 
wall and the Bogoliubov-deGennes matrix model describing the symmetry crossover CI — > 
C . The Bogoliubov-deGennes matrix model was proposed by Altland and Zirnbauer as a 
model of normal conducting-superconducting hybrid structures in mesoscopic physics [Td^ 
IT7j . It is a part of a classification scheme of random matrix ensembles in terms of the Lie 
algebra |2*2*l I23j. In the class C, the spin is conserved while the time reversal symmetry is 
broken. Then the (reduced) system Hamiltonian has a structure 

Hc=[; _" a r) (5.1) 

with an N x N hermitian a and an N x N complex symmetric b. On the other hand, 
in the class CI, the system is symmetric with respect to both spin rotations and time 
reversal. Then the Hamiltonian matrix structure is 

He, = ( I ) , (5.2) 

where a and b are both N x N real symmetric matrices. 
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Dyson proposed Brownian motion models for parametric random matrices |24j. In his 
prescription, the time evolution of a combination 



H = e- T (H C i + Ve 2 - - IK 



c 



(5.3) 



describes the crossover CI —> C. The matrix 7i is identical to Tici at r = while 7i 
approaches He as r goes to infinity. Let us redefine a and b as 



Ti 



&t -a 1 



(5.4) 



so as to reproduce (|5.1j) in the limit r — > oo. We can diagonalize 7i as 



(5.5) 



where £7 is a 2iV x 2N unitary matrix and u is a real diagonal matrix 

( ui ■ • • \ 
oj 2 '" 







\ wjv / 



(5.6) 



Assuming that 7i c is distributed according to the Gaussian distribution, we obtain the 
probability distribution function for TC as 



P{H; r)dH = A T exp 



Ti(H-e-TH C i) 



where 

N N 

dH \[ dajjdbfjdb^ Hdaf^dbf^. 

3=1 j<l 

Here the normalization constant A T is evaluated as 



(5.7) 



(5i 



(5.9) 



It is easy to calculate the differentiations of P(T~C; r) and find the Fokker-Planck equa- 



tion 



dP N 



^ KP)+ 4 (6SP)+ A fep) 



N 

+ E 

3<l 



4« p)+ 4 (4P)+ 4s (6 " p)+ 4 (6j,p) 



(5.10) 



where A is the Laplace-Beltrami operator 



A? 



A = -V 
2 frf 



5 2 



a( Qji ) 2 



«9 2 <9 2 
+ 



N 



«9 2 <9 2 
+ 



<9 2 



+ 



5 2 



<9(ajD 2 ^(4) 2 <9(^) 2 J ^ 

(5.11 
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The Laplace-Beltrami operator on a Riemannian manifold defined by the line element 

ds2 = E g^dx^dx" (5.12) 



is given by 



1 d 



det g\ dx " 



(g 1 )nv\/\det; 



I — 

dx l 



(5.13) 



In our case the corresponding line element is 

N N 

ds 2 = 2£ [(da,,) 2 + (doj) 2 + (d^-) 2 ] + 4E [(^) 2 + ( d 4) 2 + (do*) 2 + (d&j, 
j'=i j<i 

N N 

= 2^J2[da jl da* l + db jl db* l ]. (5.14) 
i=i /=i 

Let us introduce N x N matrices Ui,uz,U3 and m 4 as 

C/ = 

and rewrite da,-/ and dfyj as 



Ml «2 
M 3 M 4 



(5.15) 



N 

daji = E [(«0fci( w i)« - (usTkjius) 

k=l 
N 



Id 



dui k + dA 



■3h 



dbji 



E [{uiT kj {u 2 )ki - (U3)fej(w4)w] d^ + dB jh 

k=l 



(5.16) 



where 



N 



dAji = E [( dM i)fcj(«i)fc/ + (wi)fci(dui)w - (dM 3 )fcj(M3)fc« - (u3)fcj(dw 3 )fci 
fc=i 

d5 i« = E [( dM i)fcj( M 2)H + (ui)jy(du 2 )w - (du 3 )* kj (u A )ki - (u 3 )* kj (du 4 )ki 

k=l 

Using the unitarity of the matrix U (U*U = UW = I), we can readily see that 

N N N 

ds 2 = 2 ]T(d^) 2 + 2 ]T ]T [dA 3l dA* t + dB^ 

k=l j=l 1=1 

which, by means of ()5.12j) and (|5.13jl . yields 

. 1 * d ( T d\ . 



U k , 



UJ k . 



(5.17) 



(5.1* 



Here an operator At/ involves derivatives with respect to the variables associated with U. 



(5.19) 



The Jacobian J = J | det g | was evaluated in Ref . fUJ Ej as 

N N 

^=m^-^TiiN a - 
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Assuming that P(H; r) depends only on the radial variables u>k and a time variable r, we 
can rewrite the Fokker-Planck equation as 



Substituting P = p/ J gives 

dr ' dujj \ dojj 2 duj J 2 jr{ dujj duj 



£ = ^Ei^^^^E^Ae- (5.21) 



where 



1 * 1 

W = -J2^--logJ. (5.22) 
z i=i z 



In order to transform the Fokker-Planck operator into a form in which all the variables 
are separated, let us consider 

- ]-{H-E Q ) =e w Ce~ w (5.23) 

(Eq is a constant) and find 

N pp. N 

tf=-££l+£4 (5-24) 

For the imaginary time Schrodinger equation 

dib 1 

^ = -2^- W, (5.25) 

we call ip = G^^u^, ■ ■ ■ , wjy ; cji, • • • , Un; t) the Green function solution if it satisfies 
the initial condition 

N 

& H \uof\ ■ ■ ■ , u$>; u u ■ ■ ■ , u N ; 0) = [] ~ "f)- (5-26) 

i=i 

Since the imaginary time Schrodinger equation (J5.25j) describes the dynamics of free 
fermions, the Green function solution is given by 

G^rf', • • • , CJ!, ■ ■ ■ , u N ; r) = det[g^ (uf , u f , r)]j, l=1 ,.., N , (5.27) 

where g^ H '{u}^\ to; r) is the Green function solution of ()5.25|) with N = 1. 

Let us denote the Green function solution of the Fokker-Planck equation ()5.2()j) as 
G( FP ^ (uj^ , • • ■ , , u>i, - ■ ■ , a>jv; r). Then we can readily see that 



„ -W(ii>i,---,a;jv) 

G^)(^! 0) , • • • , wft • • • , a*; r) = e^^rG^^f, ■ ■ ■ , ■ ■ ■ , u N ; r) 

e -w(u>\ \~,uy) 



N N 

ni^r /2 e- K)2/2 n((^) 2 -M 2 ) 

e^/ 2 IT det[,((.f) 2 , (^) 2 ; r)] i)I=Wj 

i= l FT l„,(°)|i/2 p -(-< 0) ) 2 /2 TT^/.,(°)^ _ 



uy) - [u t 



n n 7 2 ri(( 

(5.28) 
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where 

OO 1 

g(x, y; r) = x^y 1 '^' 2 ^' 2 E j-C^C^y)^. (5.29) 

Here Cj(x) and hj are denned in ()3.6|) and (|3.7|) . respectively. Note that the above result 
gives the Harish-Chandra integral [23] in the case of the Bogoliubov-deGennes symmetry. 
We introduce new variables ej = up- and the corresponding Green functions as 

TV TV 

G(e?\ ..., e ( S > ;e 1 ,...,e N ;r)l[ def ) de j = G^ fp \uj[ \ • ■ ■ , w #>; ^ . . . , ^. r) n ^(0^. 

3=1 3=1 

(5.30) 

so that 

TV TV 

G(ef\ ...,«$>; Cl , .... ear; r) = e^ 2 ^ ^ de%(ef , e,; r)] jV=1 ,.. 



,7V- 



n(4 0, ) l/4 e-" /2 n( 

(5.31) 

The probability distribution functions for the eigenparameters at times r n can be 
evaluated from the Green function as 



A A . ^ . A A ■ n- ■ .AW AW. ^ \ 

V\ e li ' ' ' > e Ni T li e lJ ' ' ' ) e TV' r 2i ' ' ' i e l i'") e iV) r Mj 

de°...y o de^o(6?, e° N ) J] G(e[-\ efc 1 ; e<, e^; T(-r w ), 

(5.32) 



where po(ei, • • • , ejv) is the initial probability distribution function at r = 0. The corre- 
sponding multilevel dynamical correlation functions are given by 

n (A A ■ A A ■ ■ c M c M \ 

PBdGVtl, j fc mi i fc l! ' ' ' > t m2 ' ' ' ' ' fc l ' ' ' ' ' t m M ) 



1 ^ / C ]A ... \ AA ... A f M ... / J M 

~~ A' nAf /at \i / ue mi+l / ue TV / ue m M +l / ue TV 

^tv rn=i(^ — mi)\ Jo Jo Jo Jo 

x p(^i,---,^N;n;e 2 1 ,---,e 2 N ;r 2 ;---;ef ,■ ■ ■ ,e^;r M ). (5.33) 

Here we define the normalization constant Cat as 

poo poo poo poo 

C N = de}---/ dejv---/ def • • • / de# 

JO JO JO JO 

x P(el,---,e/V; r i; e ir--,ew;T2;---;ef ,...,ej5f;r M ). (5.34) 

The initial eigenparameter distribution for the Gaussian random matrices with the sym- 
metry ()5.2|) can be written as[T7] 

TV TV TV 

p (^i, • • • , uJN)duJi ■ ■ ■ dcjjv oc Y[ e~ auJ ^ 2 — wf \ fl |^j|dco>i ■ • • du N (5.35) 

3 3<l 3=1 



or, equivalently, 



TV TV 
-aej/2 ' 



Po(ei, • ■ ■ ,e7v)dei • • -de N oc ]\e 0163/2 ]\ Vj ~ £-l\de\ ■ ■ ■ de N . (5.36) 

3 j<i 
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A parameter a determines the variance of the Gaussian distribution for the matrix ele- 
ments and we set a = 1. 

We now see that multilevel dynamical correlation functions defined in ()5.33j) with 
the initial condition ()5.36|) have the same forms as the dynamical correlation functions for 
vicious walkers (jl.lOj) . Therefore we can similarly rewrite them in quaternion determinant 
forms 26 

PBdG(e{, ■■■,e 1 mi ;e 2 1 ,---, e 2 m2 ] • • ■ ; ef , ■ • • , e* f M ) = Tdet[5^(n /i , n v )}, 

fi,v = l,2,---,M. (5.37) 

Here the quaternion determinant is identical to that in ()3.23j) if we adopt a correspondence 

4 = A f m i = n M-i+u e 2Ti = XM-i+i- (5.38) 

Therefore all the dynamical correlation functions are shared by the matrix model and 
the vicious walk model. We have thus shown the equivalence of the vicious walk model 
with a wall in the diffusion scaling limit and the parametric Bogoliubov-deGennes matrix 
model. This equivalence holds for finite N and also in the asymptotic limits N — > oo. 
Although we have here presented the equivalence only for even N, we can similarly and 
straightforwardly prove it for odd N. 

The eigenparameter distributions of the Bogoliubov-deGennes matrix model in the 
limits r — > oo and r = are known to be equivalent to the eigenvalue distributions of 
the Laguerre unitary and orthogonal ensembles [27J |2H] of random matrices, respectively. 
Let us consider the limit r m — > oo with the time differences r m — r n fixed. We can readily 
see that in this limit the quaternion determinant is reduced to an ordinary determinant. 
The resulting determinant expressions describe temporally homogeneous dynamical cor- 
relations within the C universality class. It follows from the rescaling (J4.1)) that the C 
universality class survives until time t very close to T: only when T — t ~ 0(iV _1 ), the 
transition to CI class occurs. Therefore we can conclude that the transition from C to 
CI class is discontinuous in the limit N — > oo. The asymptotic correlation functions 
describing the CI universality class is obtained by putting v m = v n = (equivalent to 
TM-m+i = tm-u+1 = 0) in (|4.7|) . ()4.18j) and ()4.19j) . Using the Bessel function identities, 
we can easily confirm that they are identical to Nagao and Slevin's result [2Sj for the 
Laguerre orthogonal ensemble. 



6 Conclusion 

In this paper we have analyzed the vicious walk model with a wall in the diffusion scaling 
limit. It was shown that all the dynamical correlation functions are written in the forms 
of quaternion determinants. Using the quaternion determinant formulas we were able 
to derive the asymptotic formulas for the correlation functions. Finally we showed that 
the vicious walk model in the diffusion scaling limit was equivalent to the parametric 
Bogoliubov-deGennes matrix model. As the equivalence to the matrix model is so far 
established only in the diffusion scaling limit, it is interesting to consider how the matrix 
model should be generalized corresponding to the discrete vicious walk. 
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